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ABSTRACT 



Cosmic ray antiprotons provide an important probe for the study of the galac- 
tic Dark Matter, as they could be produced by neutralino annihilations, primor- 
dial black holes evaporations or other exotic sources. On the other hand, antipro- 
tons are anyway produced by standard nuclear reactions of cosmic ray nuclei on 
interstellar matter (spallations), that are known to occur in the Galaxy. This 
process is responsible for a background flux that must be carefully determined 
to estimate the detectability of an hypothetical exotic signal. 
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In this paper we provide a new evaluation of the interstellar cosmic antiproton 
flux that is fully consistent with cosmic ray nuclei in the framework of a two-zone 
diffusion model. We also study and conservatively quantify all possible sources 
of uncertainty that may affect that antiproton flux. In particular, the primary 
cosmic rays are by now so well measured that the corresponding error is removed. 
Uncertainties related to propagation are shown to range between 10% and 25%, 
depending on which part of the spectrum is considered. 



The study of the cosmic ray antiproton spectrum has been a great challenge since the 
first measurements made at the end of the seventies. Actually, the first experiments provided 
data which, in the low energy tail, showed some excess when compared to the current model 
predictions. This discrepancy stimulated a great interest into alternative explanations, viz 
the possible existence of primary antiproton sources. Such an interest did not fade even 
when further experimental data seemed to agree with theoretical predictions in standard 
leaky box models (see for example Stephens & Golden 1988 and references therein). 

Various primary antiproton sources have been proposed (Silk & Srednicki 1984; Stecker, 
Rudaz, & Walsh 1985; Ellis et al. 1988; Starkman & Vachaspati 1996; Mitsui, Maki, & Orito 
1996). The case of supersymmetric sources - relic neutralinos in the galactic halo - has 
received a particular attention and constraints on SUSY parameters have been investigated 
by comparing experimental data to theoretical predictions (Bottino et al. 1995; Chardonnet 
et al. 1996; Bottino et al. 1998; Bergstrom, Edsjo, & Ullio 1999). However, an important 
problem with this comparison is that an accurate estimation of the background secondary 
antiproton flux produced by spallations is mandatory. 

In this paper, we focus on this secondary antiproton flux, which we will call "back- 
ground" antiproton flux, having in mind the possibility of using it to determine whether 
one of the primary components ("signal") discussed above could be seen against it or not 
(P. Salati & al, in preparation; A. Barrau & al, in preparation). Such hypothetical signals 
will not be further discussed in this paper. We believe that now is a good time for a de- 
tailed evaluation of the background flux, since the next measurements of p spectra should 
be very accurate at low energy (~ 100 MeV — 10 GeV) especially in the forthcoming ten 
years (ams, bess, Pamela, . . . ). On the theoretical side, progress has already been made 
in many directions. Here are some milestones on the way: (i) the inelastic non-annihilating 
cross-section for p (Tan & Ng 1982, 1983), giving rise to the so-called tertiary contribution, 
has been taken into account (ii) the p + Hcism P contribution has been considered by 




1. 



Introduction 




3 



means of a simple geometric approach (Gaisser & Schaefer 1992), (iii) reacceleration has 
been considered (Simon & Heinbach 1996), (iv) propagation has been modeled in a more 
realistic two-zone diffusion model (Halm, Jansen, & de Niem 1993; Chardonnet et al. 1996), 
(v) the (p. He) + {H, He)isM reactions have been re-estimated in a more sophisticated nuclear 
Monte Carlo (Simon, Molnar, & Roesler 1998), (vi) the great variety of cosmic rays has been 
treated in a more coherent way (Moskalenko, Strong & Reimer 1998). As far as we know, all 
these ingredients have only been considered simultaneously in Moskalenko et al. 1998 (see 
also Moskalenko et al. 2001). 

We propose to go beyond this type of study and to use the results of our systematic 
analysis of nuclei (Maurin et al. 2001, referred hereafter as Paper I) to ascertain the theo- 
retical uncertainties on the interstellar secondary antiproton energy spectrum. This goal has 
never been achieved before, even in Moskalenko et al. 1998, 2001. The paper is organized 
as follows. Separate sections are devoted to all the ingredients entering the calculation of 
the p background: measured H and He flux, secondary production, tertiary contribution 
and propagation. Whithin each section, we first discuss the model used and the associated 
parameters; then we estimate the uncertainty they induce in the p background. An impor- 
tant aspect is worth a warning at this point. As will be discussed in section 6, the effect of 
solar modulation may be decoupled from the problem of interstellar propagation and this 
problem will not be addressed here. When a modulated flux is needed, we will use a simple 
force-field approximation modulation scheme, as in most cosmic antiprotons studies. Would 
a more careful treatment of solar modulation be needed (see for example Bieber et al. 1999), 
an interstellar flux can easily be obtained by demodulation (the force-field approximation 
modulation scheme is reversible). This interstellar flux could then be used as an input for 
any other preferred treatment of solar modulation. 

To sum up, we used results from a systematic nuclei cosmic ray analysis to consistently 
derive an antiproton secondary flux in the framework of diffusion models. As an important 
consequence we could study and quantify most of the uncertainties: in the propagation, in 
the nuclear physics and in the primary cosmic ray. We feel that our results will be valuable 
not only for speculations on primary contributions to that flux but also for the experimental 
groups which are going to perform very accurate antiproton measurements in the near future. 



2. Proton and Helium primary spectra 

The secondary antiprotons are yielded by the spallation of cosmic ray nuclei over the 
interstellar medium (see Appendix A for the formula). The most abundant species in cosmic 
rays are protons and helium, and the contribution of heavier nuclei to the antiproton pro- 
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duction is negligible. Until recently, their spectra were known with a modest accuracy and 
the data from different experiments were often incompatible at high energy. This induced an 
uncertainty of some tens of percents in the predicted antiproton spectrum. Recent measure- 
ments made by the balloon-borne spectrometer BESS (Sanuki et al. 2000) and by the AMS 
detector during the space shuttle flight (Alcaraz et al. 2000a, 2000b, 2000c) dramatically 
reduced the uncertainties both on proton and helium spectra. We fitted the high energy (T 
> 20 GeV/n) part of these measured spectra with the power law: 

$(T) = N(T/ GeV/n)-^ , (1) 

where the kinetic energy per nucleon T is given in units of GeV/n and the normalization 
factor N in units of m~^ s^^ sr"^ (GeV/n)^\ This provides a good description down to the 
threshold energy for the antiproton production. 

We fitted the BESS and AMS data both separately and combined, obtaining very similar 
results. This is obvious since the data from the two experiments are now totally compatible, 
as can be seen in Fig. 1. The upper curve presents our fit on the combined proton data. 
The best fit corresponds to = 13249 m~^ s~^ sr' ^1 (GeV/n)"' and 7 = 2.72. We do not 
plot the spectra obtained from the best fits on the single BESS and AMS data because of 
their complete overlap with the plotted curve. We did the same for helium (lower curve) and 
the corresponding numbers are = 721 m~^ s~^ sr~^ (GeV/n)~^ and 7 = 2.74. The l-cx 
deviation from the best fit spectrum does not exceed 1% for both species. Consequently, 
the corresponding uncertainty on the antiproton spectrum is smaller than the ones discussed 
in the next sections, and it will be neglected in the rest of this paper. The situation has 
significantly improved since Bottino & al. 1998, where an error of ± 25% was quoted. 



3. Antiprotons production: secondary sources 

Whereas p-p interactions are clearly the dominant process for secondary antiproton pro- 
duction in the galaxy, it has been realized long ago that p-nucleus and nucleus-nucleus col- 
lisions should also be taken into account (Gaisser & Schaefer 1992). They not only enhance 
the antiproton flux as a whole but also change its low energy tail, mostly for kinematical 
reasons. Unfortunately, very few experimental data are available on antiproton production 
cross-sections in nuclear collisions. A model-based evaluation is therefore necessary, and we 
chose to use the dtunuc program. We first discuss sub-threshold antiproton production. 
Then we present the results of our calculations of above-threshold production, which we 
compare to experimental data and analytical formulae. 
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3.1. p— p interaction 



Antiproton production via the proton-proton interaction is tlie first reaction tliat one 
has to take into account in order to evaluate the p flux. So far, the Tan & Ng parameterization 
of p cross-section (Tan & Ng 1982, 1983) has been used by almost all studies on cosmic ray 
antiprotons. To be more precise, we recall the form of secondary contribution {e.g. eq. [A9], 
Appendix A.l) 



Thus, in order to evaluate the secondary contribution of p — H^sm reaction, we used the 
parameterization of Tan & Ng (1982, 1983). We refer the interested reader to the short 
discussion in Bottino et al. (1998) for further details, or to the source papers (Tan & Ng 
1982, 1983) for a complete description. Finally, as an illustration, the impact of kinematics 
and threshold for the production rate can be found in Gaisser & Schaefer (1992). 

3.2. Calculation of the differential cross-section of antiprotons production in 



Some discrepancies between simple scalings of p-p cross-sections and experimental data 
on p-nucleus antiproton production cross-sections near threshold have been explained by 
taking into account internal nuclear Fermi motion (Shor et al. 1990). We flrst show that 
this effect does not change the cosmic antiproton spectrum. In such models, the momentum 
distribution is described by a double-gaussian function normalized to the total number of 
nucleons. The parameters are determined from scattering experiments (Moniz et al. 1971) 
and simple scaling laws. The cross-section results from a convolution 



where Pc is the internal nuclear momentum of the target nucleon, N denotes either a proton 
or a neutron (the model is isospin independent) and Ecm is the center of mass energy (with 
an off-shell target nucleon). 

Near threshold, the nucleon-nucleon cross-section can be estimated from the transition 
matrix element and the available phase space by Fermi's golden rule. Using this simple 
approach with only one free parameter (namely the matrix element), fltted on data, we have 
been able to reproduce very well most experimental results available on subthreshold antipro- 
ton production. The kinematical term was computed using a Monte Carlo multi-particle 





p-He, He-p and He-He reactions 




(3) 
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weighted event according to Lorentz-invariant Fermi phase space, whereas the integral was 
performed by adaptable gaussian quadrature. This method is not relevant to accurately 
determine the p-He, He-p or He-He cross-sections at any energy (as the momentum distri- 
bution becomes a 6 function when the involved momenta are much greater than the Fermi 
momentum) but just to investigate their behaviour below the 6 GeV kinetic energy thresh- 
old. The main result is a very fast drop below the threshold. Even after convolution with 
the ~ E'"^'^ differential power law spectrum of primary cosmic rays, two orders of magnitude 
are lost in less than two GeV below the threshold. As a consequence, the subthreshold cross 
section can be neglected to compute the secondary antiprotons flux. The above-threshold 
discrepancies between data and simple models cannot be accounted for by this effect and a 
numerical Monte Carlo approach is necessary. 

Following Simon et al. (1998), the Monte Carlo program dtunuc^ version 2.3 was 
therefore used to evaluate the cross-sections for p-He, He-p and He-He antiproton produc- 
tion reactions. The p-p reaction can be well accounted by the Tan & Ng parameterization 
(see previous section) whereas those involving nuclei heavier than helium are negligible due 
to cosmic abundances. This program is an implementation of the two-component Dual 
Parton Model (Capella et al. 1994) based on the Gribov-Glauber approach treating soft 
and hard scattering processes in a unified way. Soft processes are parameterized according 
to Regge phenomenology whereas lowest order perturbative QCD is used to simulate the 
hard component (Roesler 1997). This program calls phojet (Engel 1995) to treat individ- 
ual hadron/nucleon/photon-nucleon interaction, PYTHIA (Sjostrand 1994) for fragmentation 
of parton (according to the Lund model) and lepto ("Physics at hera" 1992) for deep 
inelastic scattering off nuclei. 



3.2.1. Comparison with experimental data 

The resulting cross-sections have been compared with experimental data on proton- 
nucleus collisions. Figure 2 shows the differential cross-section of antiprotons production 
in p-|-C and p-l-Al collisions at 12 GeV laboratory kinetic energy recently measured at the 
Proton Synchrotron in the High Energy Accelerator Research Organisation (kek-ps) for dif- 
ferent antiprotons momenta (Sugaya et al. 1998). In most cases, measurements and dtunuc 
simulations are compatible within uncertainties. The discrepancies are, anyway, taken into 
account in section 6.2 as uncertainties on the computed cross-sections. Figure 3 shows the 
invariant spectrum of antiprotons in p-|-Al collisions at 14.6 GeV/c laboratory momentum 
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as a function oi rrit — m where rrit = \/j^~-\-m^ as obtained by experiment 802 at the 
Brookhaven Tandem Alternating Gradient Synchrotron (ags) (Abbot et al. 1993). Data 
points have been normahzed by using the inelastic cross-sections and plotted for a rapidity 
interval of 1.0 < ?/ < 1.6. The results of dtunuc simulations are in perfect agreement with 
the measurements. This check is particularly important as it stands within the projectile 
energy range where most cosmic antiprotons are produced. 



3.2.2. Comparison with analytical parameterization from Mokhov & Nikitin 

Taking into account the qualitative predictions of the Regge phenomenology and partons 
model, Mokhov & Nikitin (1977) derived a parametrized inclusive cross-section for p + A ^ 
p + X: 



(Pt + ) 

where 

boPT if Pt < r ; 



' ' hj^ otherwise . 

S is the invariant mass of system, pt is the transverse momentum, xt ~ 2pT/\fS, x' = 
E*/E^^^, E* and -E^ax the total energy of the inclusive particle in the center of mass 
frame and its maximum possible value. The parameters Ci to Cq, bo, fi^ and F were not taken 
as given in Kalinovskii et al. (1989) but were re-fitted using an extensive set of experimental 
data leading to a better (Huang 2001). 

Contrary to experimental measurements that are only available for a small number 
of given energies, this analytical approach allows a useful comparison with dtunuc cross- 
sections. The resulting spectrum has therefore been propagated using the model described in 
Section 4 and the results are in excellent agreement. The dtunuc approach was nevertheless 
preferred since the Mokhov-Nikitin formula was fitted on rather heavy nuclei, and its use 
for p-He, He-p and He-He collisions would therefore require a substantial extrapolation. 



3.2.3. Results for the antiprotons 



The exclusive cross-section for antiproton production da^'^ /dEp{Ep, Ei), is obtained by 
multiplying the total inelastic cross-section of the considered reaction and the antiproton 



- 8 - 



multiplicity interaction given by dtunuc. 

This approach is time consuming since the cross-section is quite low and a large number 
of events must be generated to reach acceptable statistical uncertainties. The sampling 
points were chosen to be distributed on a logarithmic scale between 7 GeV (threshold) and 
10 TeV per nucleon for the projectile nucleus and extrapolations rely on polynomial fits. 
The antiproton kinetic energy was varied from 0.1 GeV to 100 GeV. Figure 4 gives some 
examples of differential antiproton production cross-sections as obtained from dtunuc. 



4. Tertiary contribution 

Once they have been created, antiprotons may interact with the interstellar material 
in three different ways. First, they may undergo elastic scatterings on galactic hydro- 
gen. The cross-section for that reaction has been shown to peak in the forward direction 
(Eisenhandler 1976) so that the corresponding antiproton energy loss is negligible. Antipro- 
tons are not perturbed by these elastic scatterings as they survive them while their energy 
does not change. They may also annihilate on interstellar protons. This process domi- 
nates at low energy, and its cross-section is given in Tan & Ng (1983). Last but not least, 
antiprotons may survive inelastic scatterings where the target proton is excited to a reso- 
nance. Antiprotons do not annihilate but lose a significant amount of their kinetic energy. 
Both annihilations and non-annihilating interactions contribute to the inelastic antiproton 
cross-section so that 

pp ^ pp _ pp 

^non—ann ^ me ^ ann ; V / 

where af^g is parametrized as in Tan & Ng (1983). 

For an antiproton kinetic energy Tp > 10 GeV, the Tan & Ng parameterization of 
^ann ^ which is based on experimental data - is no longer valid. The annihilation cross- 
section tends furthermore to be small at high energy. In any case, the antiproton inelastic 
but non-annihilating interaction cross section becomes equal to the total proton inelastic 
cross-section 

pp ^ pp /C-x 

^non—ann — ^ ine ' V-^/ 

The low and high energy relations for u'^on-ann do match for an antiproton kinetic energy of 
Tj5 = 13.3 GeV. 

The energy distribution of antiprotons that have undergone an inelastic but non- 
annihilating interaction has not been measured. It has been assumed here to be similar 
to the proton energy distribution after p-p inelastic scattering. An impinging antiproton 
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with kinetic energy Tp has then a differential probabihty of 

^ = A (6) 

dEp Tp 

to end up with the final energy Ep. That reaction leads to the flattening of their energy 
spectrum as the high-energy species of the peak that sits around a few GeV may replenish 
the low-energy part of the energy distribution. The corresponding source term for these 
so-called tertiary antiprotons may be expressed as 

q^i^^Ep) = / -^^{E'p^Ep}nHv'pNP{r,E'p)dE'p 

- apH^px {Ep} nnvp NP{r,Ep) . (7) 
Since the differential cross-section is given by 

the tertiary production term translates into 

f r+co pp / pi\ ^ 

gf (r,i?) = 4vrn^|y^ ^ %{r,E')dE' - <„(E) <|.,(r, E)| . (9) 

The integral over the antiproton energy E of q*p'^{E) vanishes. This mechanism does not 
actually create new antiprotons. It merely redistributes them towards lower energies and 
tends therefore to flatten their spectrum. Notice in that respect that the secondary an- 
tiproton spectrum that results from the interaction of cosmic ray protons impinging on 
interstellar helium is already fairly flat below a few GeV. Since it contributes a large fraction 
to the final result, the effect under scrutiny here may not be as large as previously thought 
(Bergstrom et al. 1999). 

As a matter of fact, antiprotons interact on both the hydrogen and helium of the Milky- 
Way ridge. Helium should also be taken into account in the discussion. As explained in 
Appendix (A. 3), we have replaced the hydrogen density in relation (9) by the geometrical 
factor riH + 4^/^ for the calculation of the tertiary component. 



5. Propagation in a diffusion model 



Propagation of cosmic rays can be studied within different theoretical frameworks, the 
most popular being the so-called Leaky Box model and the diffusion model. There is a 
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mathematical equivalence of these two approaches, which is valid only under special circum- 
stances. In particular, they lead to different results for low grammages and for unstable 
cosmic ray species (see discussion in Maurin et al. 2001). Our preference for the diffusion 
model has several justifications. First, it is a more physical approach, in the sense that 
cosmic rays are believed to diffuse in the galactic disk and halo, which is in disagreement 
with the spatial homogeneity assumed in the Leaky Box. Second, the parameters entering 
the diffusion models are related to measurable physical quantities (at least in principle), 
like the galactic magnetic field, so that their value could be cross-checked with independent 
measurements. Finally, the diffusion approach is mandatory if one wants to take primary 
sources into account, as emphasized in the introduction. 

The geometry of the problem used here is a classical cylindrical box (see for example 
Webber, Lee, & Gupta 1992) whose radial extension is i? = 20 kpc, with a disk of thickness 
2h = 200 pc and a halo of half-height L lying in the interval [1 — 15] kpc. Sources and 
interactions with matter are confined to the thin disk and diffusion which occurs throughout 
disc and halo with the same strength is independent of space coordinates. The solar system 
is located in the galactic disc {z = 0) and at a centrogalactic distance Rq = 8 kpc (Stanek & 
Garnavich 1998; Alves 2000). We emphasize that this model is exactly the one that has been 
used for the propagation of charged nuclei (Paper I) where it has been described in details. 
For the sake of completeness, we rewrite here the basic ingredients, and the parameters of 
the diffusion model we used. 



5.1. The five parameters of the model 

Our model takes into account the minimal known physical processes thought to be 
present during the propagation. Firstly, the diffusion coefficient K{E) 

K{E)=Ko(3xn^ (10) 

where the normalisation Kq is expressed in kpc^ Myr~^ and 5 is the spectral index {IZ = p/Z 
stands for the particle rigidity). Along with the spatial diffusion, one has the associated 
diffusion in energy space represented by a reacceleration term 

2 IT'S aA 

Kee{E) = -K'^=^ . (11) 
^ ^ 9 K{E) ^ ^ 

Here Kee stands for the energy diffusion coefficient which we evaluated in the no-recoil hard 
sphere scattering centers approximation. In particular Va is the alfvenic speed of scatterers 
responsible of the energetic diffusion. Next, we allow a constant convective wind directed 
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outward in the z direction. This term is represented by the velocity Vc- Motivation of such 
forms for the various parameters has been given in Paper I and will not be repeated here. 
Last, we have to include effects of energy losses. Formulae for the latter are those used for 
nuclei with the appropriated charge for an antiproton (see Paper I). 

As a consequence, diffusion model is described with five parameters: the diffusion coef- 
ficient normalization and its power index 5, the convective galactic wind velocity Vc, the 
Alfvenic speed Vj, and finally the halo thickness L. 

5.2. Configuration of the parameter space used for this analysis 

The values of these parameters are needed to compute the propagated antiproton flux. 
They may be extracted from a careful analysis of charged cosmic ray nuclei data. This has 
been done in a previous study (Paper I), where all the sets of parameters consistent with 
B/C and sub-Fe/Fe data were determined. As the propagation history for all cosmic rays 
should be similar, this is thought to be a safe procedure. In this work, we used the same sets 
and the same numerical code to propagate antiprotons, to make sure our treatment is fully 
consistent with our previous work and that the results are consistent with nuclei data. This 
is in variance with previous works using diffusion models, where the propagation parameters 
were extracted from a leaky box analysis of nuclei. It should be noticed that some of the sets 
of parameters are probably disfavored by physical considerations. For instance, our models 
have Alfven velocities Va ranging from 25 km s~^ to 85 km s~^. The upper end of this range 
is too high. Indeed, the value of the galactic magnetic field (5^1 — 2 /iG, see for example 
Han & Qiao (1994) or Rand & Lyne (1994)) and the plasma density ((^e) = 0.033 cm~^ 
according to Nordgren, Cordes & Terzian (1992)) give 10 km s~^ ^ ^ 30 km s~^. 
Besides, as mentionned in Paper I, the physical meaning of the value of Va may depend on 
the assumptions made for the scattering process. A proportionality coefficient larger than 
2/9 in relation (11) would imply smaller values for Va- The following point should also be 
kept in mind: we considered that reacceleration only occured in the thin disk, i. e. in a zone 
of half-height ha = h = 100 pc. If this process is efficient in a larger zone {ha > h), the 
overall effect is unchanged provided that the Alfven velocity is scaled down to a lower value 
as Va oc {L/haY^"^ (Seo & Ptuskin 1994). In our semi-analytical resolution of the diffusion 
model, the case ha ^ h cannot be straightforwardly taken into account, but the previous 
conclusion would still hold. Indeed, we can make the reacceleration zone larger by increasing 
the disk thickness h, while keeping constant the quantity nnh so that all the other effects are 
unaffected. For example, a ha = 1 kpc reacceleration zone would lead to Alfven velocities 
about three times smaller so that in the sets of parameters used in this study, Va would range 
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between ~ 10 km s ^ and ~ 30 km s ^. Anyway, we adopt a conservative attitude and we 
do not apply any cut in our initial sets of parameters. 

To sum up, we have applied all the configurations giving a good (less than 40 for 26 
data points and 5 parameters) in the B/C analysis of Paper I (see this paper for an extensive 
description of the nuclei analysis). We insist on the fact that none of this parameter is 
further modified or adjusted, they are not free parameters. 

5.3. Calculation of the secondary component 

Once the set of diffusion-propagation parameters is chosen as explained above, evalua- 
tion of the corresponding flux is straightforward. Semi-analytical solution for the antiproton 
background is given in Appendix A. Apart from the propagation, the two other necessary 
inputs are - as one can see from equation (A9) - the measured top of atmosphere H and He 
flux discussed in Section 2, and the nuclear processes described in Section 3 and 4. 

To compare our results to experimental data, solar modulation (the effect of the solar 
wind on the interstellar flux crossing the heliosphere) must be taken into account. We chose 
to use the so-called force-fleld approximation, which is used in most antiproton studies (see 
last section for a discussion). 

In all the subsequent results, the top-of-atmosphere antiproton flux has been obtained 
from the interstellar one with a modulation parameter of = 250 MV ($ = Z/A x = 
250 MV), adapted for a period of minimal solar activity. This choice is motivated by the 
comparison to BESS data taken during the last solar minimum. 

6. Results and uncertainties 
6.1. Results 

We have calculated the secondary, top-of-atmosphere antiproton spectrum obtained 
with the procedure described above. To begin with, we chose a particular set of diffusion 
parameters giving a good flt to the B/C data (see above). Namely, we have flxed: Kq/ L = 
0.0345 kpc^ Myr~^, L = 9.5 kpc, Vc = 10.5 km/s and Va = 85.1 km/s. This set gives the best 
for 6 flxed to 0.6 and the resulting antiproton spectrum will be used as a reference in most 
subsequent flgures. Fig. 5, displays this computed antiproton flux along with experimental 
data collected by the BESS spectrometer during two flights in a period of minimal solar 
activity. Circles correspond to the combined 1995 and 1997 data (Orito et al. 2000) and 
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squares to the 1998 ones (Maeno et al. 2000). The dotted hnes represent the contribution to 
the total flux coming from the various nuclear reactions: from top to bottom are represented 
the contribution of p-p, p-He, He-p and He-He. 

First of all, we notice that the calculated spectrum agrees very well with the BESS data 
points. This strong result gives confidence in our consistent treatment of nuclei and antipro- 
ton propagation. Second, even if the main production channel is the spallation of cosmic 
ray protons over interstellar hydrogen, we see that the contribution of protons over helium is 
very important, particularly at low energies (where a hypothetical primary signature would 
be expected). It emphasizes the necessity of having a good parameterization of the p-He 
reaction. 

In the following sections, we study and quantify all the uncertainties and possible sources 
of errors in the secondary antiproton flux given above. 

6.2. Uncertainties from diffusion parameters 

The first source of uncertainty comes from the fact that the propagation parameters are 
not perfectly known, even if they are severely constrained by the analysis of B/C experimental 
results (Paper I). A quantitative estimate for this uncertainty is obtained by applying all the 
good parameter sets to antiproton propagation. In a first step, we set the diffusion coefficient 
spectral index 6 to 0.6 and allow the four other parameters {Kq, L, Vc and Va) to vary in 
the part of the parameter space giving a good fit to B/C. The resulting antiproton fiuxes 
are presented in Fig. 6. The two curves represent the minimal and the maximal fiux obtain 
with this set of parameters. In a second step, we also let 6 vary in the allowed region of 
the parameter space, along with the four other parameters (Fig.7 and Fig.8 of Paper I). As 
before, the minimal and maximal fiuxes are displayed in Fig. 7. The resulting scatter depends 
on the energy. More precisely, it is 9% from 100 MeV to 1 GeV, reaches a maximum of 24% 
at 10 GeV and decreases to 10% at 100 GeV. This gives our estimate of the uncertainties 
related to diffusion. They may be considered as quite conservative, as the range of allowed 
parameters could probably be further reduced by a thorough analysis of radioactive nuclei 
(Donato et al, in preparation) and also by new measurements of stable species. 

6.3. Uncertainties from nuclear parameters 

The uncertainties on the antiproton production cross-sections from p-He, He-p and 
He-He reactions have been evaluated using the most extensive set of experimental data 
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available. In addition to those described in Section 3.2.1, the average antiproton multiphcity 
in p-p coUisions as measured by Antinucci et aL (1973) has also been checked out. Finally, 
measurements from Eichten et al. (1972) performed by the CERN-Rome group with the 
single-arm magnetic spectrometer (AUaby et al. 1971) were taken into account. They give 
the Lorentz invariant density (defined as 2ES'^a / (aap'^SpSQ) where E and p are the laboratory 
energy and momentum of the produced antiproton and cJa is the absorption cross-section) 
as a function of p and of the production angle 6. A wide range of values from 6 = 17 mrad 
to ^ = 127 mrad and from p = A GeV to p = 16 GeV has been explored. 

All those measurements have been compared with dtunuc results. As mentionned 
before, most of them are in excellent agreement with the simulation. The more important 
discrepancies were found for high-energy produced antiprotons in p-Be collisions and for low 
energy projectile protons in p-p collisions. This latter point is not surprising as the physical 
input of DTUNUC can hardly be justified for a center of mass energy ^/s < 10 GeV. In 
both cases, experimental cross-sections were lower than the simulated ones. Differences are 
never larger than a factor of two. To account for such effects we parameterized maxima and 
mimina cross-sections as a correction to the computed ones, depending on the projectile and 
antiproton energies. The simplest, i.e. linear, energy variation was assumed and the slope 
was chosen to be very conservative with respect to experimental data. Finally, it has been 
checked that changes in the Monte Carlo results induced by small variations of the input 
physical parameters remain within the previously computed errors. 

According to Tan & Ng (1982, 1983), the uncertainty in the parameterization of their 
p-p cross-section should not exceed 10%. From another point of view, Simon et al. (1998) 
have compared two parameterization of the existing data along with the Monte Carlo model 
DTUNUC. They found large discrepancies which induce a 40% effect on the antiproton 
prediction. Nevertheless, since data are available for that reaction, we think that the Tan & 
Ng parameterization is more reliable than any Monte Carlo. 

In Fig. 8 we present our estimation of the uncertainties related to nuclear physics. 
The central curve is our reference presented above. The upper one is obtained with the 
set of maximal p-He, He-p, He-He cross-sections while increasing the p-p cross-section by 
10%. Similarly, the lower curve is obtained with the minimal values for these cross-sections 
while decreasing the p-p cross-section by 10%. Indeed, such a variation for p-p has been 
included for the sake of completeness even if it modifies the antiproton spectrum only by a 
few percents. As a conclusion, the shift of the upper and the lower curve with respect to the 
central one is of the order of 22-25 % over the energy range 0.1-100 GeV. 

Besides these major sources of uncertainties, we have also investigated the influence 
of a possible error in the parameterization of the inelastic non-annihilating cross-section. 
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which gives rise to the tertiary component. We modified it by 20%, which is thought to 
be very conservative. We found that the antiproton spectrum is modified by less than 1%. 
In the same fine of thought, the effect of total inelastic plus non-annihilating reactions on 
interstellar He is found to be negligible (see discussion in Appendix A. 3). 

6.4. Other uncertainties 

There are few other sources of uncertainties. To begin with, as we discussed in Section 
2, primary cosmic ray fluxes (protons and helium) have been measured with unprecedented 
accuracy. For the flrst time, the induced uncertainties on the antiproton spectrum can be 
neglected. 

Next, the only parameters which have not been varied in the previous discussion are 
those related to the description of the interstellar medium, i.e. the densities nn and riHe- In 
all the preceding analysis, these were flxed to n^su = "nn+n-nc = 1 cm~^ and fuc = nuc/nisM = 
10% (same as in Paper I). We have tested the sensitivity of our results to changes in both 
riisM and /hc- For this purpose, we found the new values for the diffusion parameters (for 
6 = 0.6) giving a good flt to B/C, and applied them to antiprotons. Varying /hs in the range 
5% < /hc < 15%, the resulting flux is modifled by less than 15% over the whole energy range. 
Notice that this range of /hc values can be considered as very conservative (see discussion 
in Strong & Moskalenko 1998). A more reahstic 10% error on fue (i-e. 0.9% < fuc < 1-1%) 
would lead to a few % error on the antiproton spectrum. Alternatively, varying njsM from 
0.8 to 1.2 cm~^, the resulting flux is modifled by less than 0.5% over the whole energy 
range. To sum up, the only contributing errors are from the helium fraction /hc through the 
dependence of antiproton production on corresponding cross-sections. 

Finally, solar modulation induces some uncertainty. This problem is still debated, and a 
rigorous treatment of this effect is beyond the scope of this paper (see for example Bieber et 
al. 1999 for a recent analysis). However, in a "force-fleld" approximation, a general feature 
is that the steeper the spectrum, the greater the effect. Our antiproton spectra being rather 
flat, we do not expect them to be dramatically affected by a change in the modulation 
parameter. Anyway, this local effect is decorrelated from the propagation history. Solar 
modulation - which is the last energetic modiflcation suffered by an incoming galactic cosmic 
ray - can thus be treated completely independently from the above analysis. Figure 9 shows 
our demodulated spectra together with other interstellar published spectra (Simon et al. 
1998, Bieber et al. 1999, Moskalenko et al. 2001) 



- 16 - 



7. Conclusions 

We have computed cosmic antiproton fluxes in the framework of a two-zone diffusion 
model taking into account galactic wind, stochastic reacceleration and energy losses. The 
propagation parameters have been chosen according to Maurin et al. (2001), as to be in agree- 
ment with cosmic ray nuclei data. The annihilating as well as the inelastic non-annihilating 
(tertiary) p-p reactions have been taken into account. The p-p, He-p, p-He and He-He 
nuclear reaction have also been included and the relevant cross-sections have been computed 
using the Monte Carlo program dtunuc. The latest measured values for cosmic protons 
and helium fluxes from AMS and BESS have been considered. 

The results may be summarized as follows. First, the values of all the inputs being 
either extracted from the analysis of nuclei (diffusion parameters 6, L, Kq, Vc and Va) or 
measured (proton and helium fluxes), all the cosmic antiproton fluxes naturally coming out 
of the calculation are completely contained within the experimental error bars of BESS data. 

The other strong conclusion is that all possible sources of uncertainties have been de- 
rived. They have been signiflcantly improved with respect to the previous gross estimates. 

In particular, those related to propagation range between 10% and 25%, depending on 
which part of the spectrum is considered, and those related to nuclear physics are below 
25 %. We emphasize that the uncertainties related to propagation will probably be further 
reduced by a more complete study of cosmic ray nuclei, in particular by focusing on the 
radioactive species. We also note that more accurate data on cosmic ray nuclei fluxes would 
give better constraints on the diffusion parameters, which in turn would translate into lower 
uncertainties on antiprotons fluxes. The major remaining uncertainties come from nuclear 
physics and are already comparable to experimental error bars. As antiproton spectrum 
measurements should better in the near future, antiproton studies could be limited by nuclear 
undeterminacies. Further work and especially new measurements of antiproton production 
in the p-He channel would be of great interest. 
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A. Solution for the secondary antiprotons 

We summarize in this annex the salient features of our derivation of the spallation 
antiproton energy spectrum. The propagation of cosmic-rays throughout the galaxy is de- 
scribed with a two-zone effective diffusion model which has been thoroughly discussed in 
a preceeding analysis (Paper I). The Milky- Way is pictured as a thin gaseous disk with 
radius R = 20 kpc and thickness 2h = 200 pc where charged nuclei are accelerated and 
scatter on the interstellar gas to produce in particular secondary antiprotons. That thin 
ridge is sandwiched by two thick confinement layers. The effective diffusion of cosmic-rays 
throughout the galactic magnetic fields occurs uniformly within the disk and halo with the 
same strength. Furthermore, we consider here a constant wind Vc in the z direction. The 
associated adiabatic losses take place in the disk only. 



A.l. High energy limit 

As compared to the cosmic-ray nuclei on which the analysis of Paper I has focused, 
antiprotons have the same propagation history but differ as regards their production. The 
space-energy density A^^ is related to the antiproton flux through 

$p-(r,i?) = -^Vp{E)NP{r,E) . (Al) 

As explained in Paper I - see in particular their equation (Al) - the density A^^ satisfies the 
relation 

2hS{z)q'/'{r,0,E) = 2hSiz)Tf'NP{r,0,E) + 

as long as steady state holds. Diffusion and convection have been included. Inelastic in- 
teractions on interstellar atoms are described through the collision rate P^"*^ which will be 
discussed in more detail together with tertiary antiprotons. The antiproton density 

CXD 

N^{r,z,E) = Y,N!{z,E)Jo{c~) , (A3) 

i=l 

and the secondary source term 

oo 

g|-(r,0,E) = J2<lp7iE)Jo{Q^) , (A4) 
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may be expanded over the orthogonal set of Bessel functions JqIQx) where Q stands for the 
ith zero of Jq while z = 1 . . . oo. The boundary condition A^^ = is therefore readily ensured 
for r = R. The Bessel transform of the antiproton density has vertical dependence 

Nf{z,E) = Nf{E) expj^l {sinh || (L - .) | / sinh || l| | , (A5) 



where the quantity Si is defined as 



ry2 ^2 >| 1/2 



Solving equation (A2) with the help of the Bessel expansions (A3) and (A4) leads to the 
simple relation 

Nf{E) = ^ , (A7) 

that mostly holds at high energy - say above ~ 100 GeV - where energy losses and diffusive 
reacceleration do not play any major role. The coefficients are given by 



Af{E) = 2 h Vp"' + V, + K S, coth . (A8) 

Notice that the diffusion coefficient K - that comes into play in the definition of Si and 
therefore of - essentially depends on the rigidity. One should keep in mind that the 
relationship between K and the energy per nucleon may actually depend on the nuclear 
species at stake through the average charge per nucleon Z/A. 

Secondary antiprotons are produced by the spallation reactions of high-energy cosmic-ray 
protons and helium on the interstellar material of the Milky Way ridge at 2; = 0. The source 
term 

ql'^r.E) = / —{a{E') + h^p{E)}n,{An^,{r,E')}dE' (A9) 

^Threshold "-^^ 

corresponds to particles a - protons or helium - impinging on atoms h - hydrogen or helium 
- at rest. Four different production channels need therefore to be considered depending on 
the nature of the cosmic-rays and of the stellar gas. Proton-proton collisions are discussed 
in Section (3.1) whereas interactions that involve at least a helium nucleus are reviewed in 
Section (3.2). Bessel expanding relation (A9) leads to 

poo J 

qpm= / —{a{E') + b^p{E)}n,v'NnE')dE'. (AlO) 

J Threshold "-^^ 

The primary species a are accelerated in the galactic disk so that their own production rate 
may be expressed as 

(r, z,E) = 2h S{z) q, (r, 0, E) oc (i?) /(r) , (All) 
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where q^°^{E) denotes the global galactic production rate of particles - protons or helium - 
with energy E in the energy bin dE. We have assumed here that the energy dependence 
of that production rate could be disentangled from its distribution /(r) along the galactic 
disk. The bulk of the secondary antiproton production takes place for a typical energy of 
the impinging species of ~ 20 — 30 GeV/n. Note also that the primary fluxes ^a{E) are 
monotonically decreasing with the energy E. Both energy losses and diffusive reacceleration 
have therefore a negligible effect on the spectra $„. We readily infer that the Bessel transform 
A''" may be expressed as 

= C{E) , (A12) 

where the coefficients A"^ are given by a relation similar to (A8) whereas the quantities qi 
are defined as 

" jfTCO {/ '^duJoiCiu} f {r = uR}^ l^j^ ud%if{r = uR}^ . (A13) 

The cosmic-ray flux may be determined everywhere as it is related to the Bessel transform 
N"" through relations similar to (Al) and (A3). The cosmic-ray flux scales in particular 
with the global galactic production rate This allows to determine the latter by imposing 
that the interstellar proton and helium fluxes at the solar system do actually match the 
observations. 



A. 2. Full solution without tertiaries 

Forgetting for a while that the inelastic collisions of antiprotons with the interstellar 
gas may be disentangled into annihilating and non-annihilating interactions, we have to 
modify relation (A2) so as to take into account now the energy losses as well as diffusive 
reacceleration. This is straightforward since those processes take place only in the disk and 
not in the halo. Once again, following the procedure described in Paper I, one gets the 
differential equation 

ATf + 2hdE{hUE)Nf - KUE)dENf} = 2hq;r{E) , (A14) 

where fe^^^^ and K^^ stand respectively for the energy losses and the diffusion in energy. 

A. 3. Full solution with tertiaries 

We have seen that the source term for tertiaries is 
qf{r,E) = 47rnH|y^ ""^ ' ^,{r,E')dE' - $,(r, E)| . (A15) 
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Remembering that the antiproton flux $p is related to the space-energy density A^^ through 
equation (Al) and Bessel expanding relation (A15) leads to 

qf-iE) = ' nHv'Nf{E')dE' - aZn-ann{E) uh v Nf{E) . (A16) 

In the thin disk approximation, that expression needs to be multiplied by 2h6{z). The 
Bessel transforms Nf{z = 0,E) of the antiproton density obey now the integro-differential 
equation 

A^Nf + 2hdE{hLs{E)N! - Ki^{E)dENf] = 2h {ql7{E) + qfKE)} • (A17) 

Notice that in the deflnition of the coefficients , the rate r*^"^ should now be replaced by 

(E) = aZniE)v,{E)nH , (A18) 

where annihilations alone are considered. The inelastic non-annihilating reactions are di- 
rectly dealt with in the tertiary production term qp^. 

Helium should also be taken into account in our discussion of the annihilations as well as 
of the inelastic but non-annihilating interactions which antiprotons undergo with interstellar 
matter. As there are no measurements, we have adopted as an educated guess the geometrical 
approximation which consists in scaling the appropriate cross-sections by a factor of 4^/'^ when 
we deal with helium. In the formulae (A16), (A16) and (A18), we have therefore replaced the 
hydrogen density uh by (uh + 4^/^ nne). Such a replacement has httle effect. That overall 
change in the propagated antiproton spectrum is at most 1%. 



B. Numerical resolution 

We need now to solve the energy-diffusion equation (A17) for each Bessel order i. In the 
absence of diffusive reacceleration and energy losses, its solution Nf ° satisfies the relation 

A^Nf' = 2h{qir{E) + qfKE)} ■ (Bl) 

Defining the functions 

CiE) = , (B2) 

and 

a{E) = , (B3) 
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where T = Ep — nip is the antiproton kinetic energy, allows us to simplify equation (A17) 
into 

« + c4-lbLu - = , (B4) 



where and u respectively stand for Nf^ and Nf . We can express relation (B4) on a 
one-dimensional grid extending from Xjnf to Xsup with x = In (T/Tjnf). We are interested in 
kinetic energies extending from Tinf = 100 MeV up to Tsup = 100 GeV. The spacing between 
two points in energy is 

Ax = ^ In , (B5) 



N T, 



' mm 



where N has been fixed to 150 in our code. Our resolution method lies on the direct inversion 
of the algebraic linear equations that translate relation (B4) on the set of the + 1 different 
values of the variable x. If j denotes the point at position 

Xj 

we get 



The matrix A that connects u to has been written here so as to be tridiagonal. This 
allows for a fast inversion of the algebraic equation (B7). 

• For < j < N, the tridiagonal matrix A may be written as 



' 2Ax^'-' Ax 

while 



C C 

^jj-i = -irj-i>"-i - , (B8) 



bj = = 1 + {aj-i/2 + ai+1/2) , (B9) 



and 



C C 

= = J^bf^i - , (BIO) 



• The boundary j = corresponds to the low energy tip Tmin = 100 MeV where we have 
implemented the condition u(xmin) = 0. This translates into U-1/2 = ^1/2 and leads to the 
matrix elements 

ao = Ao,-i = , (Bll) 

and 

60 = Ao,o = 1 - |^^o°" + ^ («l/2-«-l/2) , (B12) 



- 22 - 



and also 



Co = Ao,i = - ^ (ai/2-a_i/2) . (B13) 



• We have finally assumed that both u and m° were equal at the high-energy boundary 
j = N. In this regime, the energy losses and the diffusive reacceleration should not affect 
too much the cosmic-ray energy spectrum. This translates into the simple conditions 

aN = An,n-i = , (B14) 

and 

bN = An,n = 1 , (B15) 

whereas, by definition 

Cat = ^AT.TV+l = . (B16) 

Inverting a tridiagonal matrix such as A may be potentially dangerous as Jordan pivoting 
is not implemented in the standard resolution scheme. As a matter of fact, energy losses 
and diffusive reacceleration lead to a moderate change in the antiproton spectrum. This 
translates into the fact that the matrix A is close to unity. We have nevertheless checked 
that our results remained unchanged when Gauss- Jordan inversion was used (Press et al. 
1992). We have also modified relation (B4) into the time-dependent equation 



a^ + " + C- 6,L"-«;^^ =0. (B17) 



It may be shown that the static solution u to equation (B4) also obtains from the superpo- 
sition 



-oo 

burst/' 



u = / u'"'''\t)dt , (B18) 
Jo 

of the reaction to an initial burst 

yburst^^-^ = n° (B19) 

taking place at t = and subsequently evolving according to relation (BIT). The later 
equation has also been solved on a discrete set of + 1 values of the antiproton kinetic 
energy while a Crank-Nicholson scheme was implemented. Once again, the result (B18) is 
the same as what the direct inversion of the algebraic set of relations (B7) gives. We are 
therefore confident that our resolution procedure is robust. 



The tertiary source term depends on the global antiproton energy spectrum that is itself 
determined by the differential equation (A17). Starting from a trial antiproton spectrum - 
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say for instance Nf^ with only the secondary production mechanism qp^^[E) cranked up - we 
invert equation (B7). The new energy spectrum is used to compute the tertiary source term 
q*pl{E) through the integral (A16). We may therefore proceed once again through the same 
steps and invert the diffusive reacceleration equation (A17) until the antiproton spectrum 
becomes stable. We have actually checked that convergence obtains after ~ 5 recursions. 
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Fig. 1. — The upper (lower) curve displays the measured proton (helium) flux along with an 
analytical flt (see text). On both curves, data are from AMS (Alcaraz et al. 2000a, 2000b, 
2000c) (crosses) and BESS (Sanuki et al. 2000) (fllled circles). 
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Fig. 2. — Here are displayed the antiproton production cross-section in p+C (top) and p+Al 
(bottom) collisions at 12 GeV laboratory kinetic energy. Filled circles are experimental data 
(Sugaya et al. 1998) and the lines are from our dtunuc simulations. The error bars have 
been assumed to be 15%. This value is usual for such experiments and was suggested by a 
analysis combining most data available. 
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Fig. 3. — Invariant spectrum of p in p+Al collisions at 14.6 GeV laboratory momentum. 
Filled circles are experimental data (Abbott et al. 1993) and the line is from our dtunuc 
simulation. 
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Fig. 4. — From top to bottom: antiproton differential production cross-section in He-He, 
p-He and He-p reactions for antiprotons kinetic energy 1.5 GeV, as obtained witli dtunuc 
simulations. 
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Fig. 5. — Solid line shows the total top-of-atmosphere (toa) secondary antiproton spectrum 
for the reference set of diffusion parameters (see text for details). Dashed lines are the 
contributions to this total flux from various nuclear reactions (from top to bottom: p-p, 
p-He, He-p and He-He). Data points are taken from BESS 95+97 (filled circles) and from 
BESS 98 (empty squares). 
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Fig. 6. — This plot shows the envelope of the TOA antiproton spectra generated with the 
sets of diffusion parameters consistent with B/C and for which 6 has been fixed to 0.6 (data 
points are the same as in Fig. 5). 
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Fig. 7. — Same as previous figure, but where the whole region of parameter space consistent 
with B/C has been used (Fig. 7 of Paper I). The resulting bounds give an estimation of the 
uncertainty due to the undeterminacy of the diffusion parameters (data are the same as in 
Fig. 5). 
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Fig. 8. — In this figure, the toa antiproton spectrum has been computed with extreme 
values of dtunuc nuclear parameters. The central line is the reference curve showed in Fig. 
4, while upper and lower curves correspond respectively to the maximum and minimum of 
the antiproton production rate. These two bounds give an estimation of the uncertainty due 
to the undeterminacy of the nuclear parameters (data are the same as in Fig. 5). 
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Fig. 9. — Comparison of our interstellar spectra (thick solid lines indicate the lower and 
upper band due to the uncertainties in the propagation parameters) with other published 
antiproton spectra. Dotted lines are lower and upper values from Simon et al. (1998), short 
dashed line is from Bieber et al. (1999) and long dashed line is from Moskalenko et al. 
(2001). 



